Gutzwiller description of non-magnetic Mott insulators: a dimer lattice model. 
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We introduce a novel extension of the Gutzwiller variational wavefunction able to deal with 
insulators that escape any mean-field like description, as for instance non-magnetic insulators. As 
an application, we study the Mott transition from a paramagnetic metal into a non-magnetic Peierls, 
or valence-bond, Mott insulator. We analyze this model by means of our Gutzwiller wavefunction 
analytically in the limit of large coordination lattices, where we find that: (1) the Mott transition 
is first order; (2) the Peierls gap is large in the Mott insulator, although it is mainly contributed by 
the electron repulsion; (3) singlet-superconductivity arises around the transition. 

PACS numbers: 71.10.-w, 71.10.Fd, 71.30.+h 



I. INTRODUCTION 

Among the theoretical tools devised to deal with 
strongly correlated metals close to a Mott metal-to- 
insulator transition (MIT), the simplest one likely is 
the variational approach introduced in the 60ths' by 
Gutzwiller— £ to describe itinerant ferromagnetism and 
narrow band conductors. In its original version, the 
Gutzwiller variational wavefunction has the form 



(i) 



R 



where \<p) is an uncorrelated wavefunction for which 
Wick's theorem holds, Vr. an operator at site R, and 
both |0) and Vr, have to be determined by minimizing 
the variational energy. The role of the operator V-r. is to 
modify, according to the on-site interaction, the weights 
of the local electronic configurations with respect to their 
values in the uncorrelated wavefunction. 

In spite of its simplicity, the Gutzwiller wavefunction 
is quite effective in capturing physical properties that 
supposedly identify strongly correlated metals, as for in- 
stance the large increase of the effective mass. 3 However, 
since the dependence upon the distance |R — R'| of inter- 
site correlations are still determined by the uncorrelated 
wavefunction, while the local operators Vr just affect the 
amplitudes, the Gutzwiller wavefunction can describe a 
Mott insulator either if Vr suppresses completely charge 
fluctuations, that provides a very poor description of an 
insulator, or if \(f>) itself is insulating. The latter case 
can be stabilized within the original Gutzwiller approach 
only when \<f>) is an admissible Hartree-Fock solution of 
the Hamiltonian. As an example let us consider a sin- 
gle band model at half-filling, for instance the Hubbard 
model 



H 



RR'.( 



£rr' 4 ff c R , 



n R)T n R ,j., 



r 



where c R a and c R a creates and annihilates, respectively, 
an electron with spin a =T,I at site R and riR^ a = 



ct, Cx, . This Hamiltonian admits at the mean-field 
level two possible phases, one paramagnetic, (tir^) — 
(nR,j), and the other magnetic, (nRj) ^ ( n R,l)- The 
latter is the only one that can eventually describe an 
insulator. In the Hubbard model the action of the oper- 
ator V is to increase the weight of singly occupied sites 
at expenses of doubly occupied and empty sites, in or- 
der to minimize the Coulomb repulsion U. Evidently, 
even when the repulsion is very strong, hence the model 
is a Mott insulator, a realistic wavefunction should still 
allow for charge fluctuations responsible for the super- 
exchange, that survives even deep inside the Mott phase. 
However, any paramagnetic uncorrelated wavefunction, 
for instance the Fermi sea, is unable to generate any 
super-exchange and necessarily leads to a non-realistic 
Mott insulator where configurations with empty or dou- 
bly occupied sites are fully suppressed.— The only way 
to generate super-exchange is to assume a magnetically 
ordered \<p), which is also the only insulating wavefunc- 
tion accessible within Hartree-Fock. However a magnetic 
state might not always be the right choice, especially if 
magnetism is sufficiently frustrated. 

Recently, an improved version of the Gutzwiller wave- 
function has been proposed,— in which additional inter- 
site correlations are provided by density-density Jastrow 
factors, namely 



G ) -> exp [ - ^ W R,R' n R n R' 



RR' 



(2) 



where tir is the site R occupation number and wr ; r' 
variational parameters. This novel class of wavefunctions 
has the capability to disentangle charge from other de- 
grees of freedom, hence is more suitable to capture Mott 
localization, as it has indeed been shown^^ However, 
unlike the conventional Gutzwiller wavefunction (JTJ) , the 
Gutzwiller- Jastrow wavefunction ^ can only be dealt 
with numerically by variational Monte Carlo, which is in- 
herently limited to finite-size systems, albeit quite large^ 
An alternative approach, that is closely related to 
recently proposed extensions of Dynamical Mean Field 
Theory (DMFT) from the original single-site formula- 
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tion^ to a cluster on o 9 ' 10 i n ' 12 i 13 ; is to consider a vari- 
ational wavefunction of the same form as ([T|) but de- 
fined on a lattice with non-primitive unit cells. In this 
case, the operator Pr acts on all the available electronic 
configurations of the lattice sites belonging to the non- 
primitive cell. The advantage is that in this way one may 
include additional short-range correlations without losing 
the property of the wavefunction to be analytically man- 
ageable, at least in infinite-coordination lattices. The ob- 
vious disadvantage is that this wavefunction could bias 
the variational solution towards translational-symmetry 
breaking. 

Within this scheme, the variational problem becomes 
generically equivalent to optimize a Gutzwiller wave- 
function for a multi-band Hamiltonian. There have 
been recently an amount of attempts to extend the 
Gutzwiller wavefunction to multi-orbital models that in- 
clude further complications like for instance Coulomb ex- 
chang o 14 ! 15 ' 16 ! 17 . In this paper we introduce a further 
extension that is capable to generate inter-site correla- 
tions as the super-exchange for paramagnetic wavefunc- 
tions, otherwise missed by the conventional Gutzwiller 
approach. This novel class of wavefunctions also allows 
to explore new kinds of variational solutions. Specifi- 
cally, there are interesting examples of correlated models 
where the Mott insulating phase escapes any Hartree- 
Fock mean-field treatment, in other words can not be 
represented by a single Slater determinant. A very sim- 
ple case, that we will explicitly consider throughout this 
work, is a Peierls insulator, namely a short-range valence- 
bond crystal, in which pairs of nearest neighbor sites are 
strongly bound into a singlet configuration, leading to a 
state that is simply a collection of spin-singlets. Such a 
Mott insulating state is not accessible by Hartree-Fock 
theory, just because each singlet is itself not express- 
ible as a Slater determinant, nor by the conventional 
Gutzwiller approach, which, as mentioned, gives a poor 
description of paramagnetic insulators. 

The paper is organized as follows. In Section |TT] we 
present the variational wavefunction and discuss under 
which conditions it can be deal with analytically. In Sec- 
tion IIIII we discuss how to build up the wavefunction in 
the case in which the basic unit of the lattice model is 
a dimer. Next, in Section IIV1 we solve the variational 
problem for a specific lattice model of dimers. Conclu- 
sions are given in Section fVl 



II. THE VARIATIONAL WAVEFUNCTION 

In this Section, we introduce an extension of the 
Gutzwiller wavefunction (JTJ) which is particularly con- 
venient to perform analytical calculation in the limit of 
infinite-coordination lattices.— Let us consider a generic 
multi-band Hamiltonian. Each lattice site R contains 
several orbitals that give rise to a bunch of electronic con- 
figurations which we denote individually as |r;R). The 



most general operator Pr can be chosen of the form: 

Vn = J2 A(R)rr'|r;R)<r';R|, (3) 

IT' 

where A(R)rr' are variational parameters. In general 
needs not to be hermitean, namely for r ^ T' it is not 
required that A(R) rr , = A(R)r'r- Indeed, as we shall 
see, the non-hermitean character plays a very important 
role. We further assume that the Wick's theorem holds 
for the uncorrelated wavefunction, hence that \<f>) is either 
a Slater determinant or a BCS wavefunction. 

It was realized by Biinemann, Weber and Gebhard 1 ^ 
that average values of operators on the Gutzwiller wave- 
function (jTJ) can be analytically computed in infinite co- 
ordination lattices provided the following two constraints 
are imposed on Pr: 

M^r» = <^> =1. ( 4 ) 
(0I^r^rCr|0> = MC R |0>, (5) 

where Cr is the local single-particle density-matrix op- 
erator, with elements c R a c R p and c R a c R p, a labeling 

single-particle states, while c R and c R a create and an- 
nihilate, respectively, an electron at site R in state a. 

The first constraint, Eq. (T5|), does not actually limit 
the variational freedom, since Pr is defined up to a nor- 
malization factor. On the contrary, the latter constraint, 
Eq. ([5]), may reduce the variational freedom, although it 
seems not in a relevant manner, at least in all cases that 
we have so far investigated. We notice that Eq. is not 
the same as imposing 

(0|PrCrP r |0) = (0|Cr|^ (6) 

unless P R commutes with Cr, which is a further con- 
straint to be imposed on Pr. This actually is the only 
case that has been hitherto considered, see e.g. Refs. \l4 
and [H. However, as we shall see, there are interesting 
models which force to abandon the supplementary con- 
dition ([6]), which is anyway unnecessary^ 

By means of Wick's theorem, the left-hand side of © 
includes a disconnected term 

MPrPh» W>|Cr|0) = MCrM, 

where the right-hand side follows from Q, plus con- 
nected terms that are obtained by selecting in all possible 
ways a pair of single- fcrmion operators from P R P R , av- 
eraging on \(j)) what remains, and finally averaging the 
two single-fermion operators with those of Cr. There- 
fore, imposing (T5J) means that the sum of all connected 
terms vanishes, whatever is the element of the single- 
particle density-matrix. In other words, the operator 
that is left after taking out from P R P R any pair of single- 
fermion operators has null average on \<p). In turns, this 
also implies that, when averaging on \4>) P R P R with 
multi-particle operators at different sites, the only con- 
nected terms that survive are those that involve four or 
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FIG. 1: Graphical representation of the average on \cj>) of 
V^Vji, drawn as a box, times a generic multi-particle oper- 
ator at site R', drawn as a circle. Lines that join the two 
operators represent the average of two single-fermion opera- 
tors, one at R and the other at R'. The dots include all terms 
where the two sites are joined by more than four lines. The 
important thing to notice is the absence of terms in which the 
two sites are connected by two lines. 

more single-fermion operators of P R P R , that are repre- 
sented graphically in Fig. [1] as lines coming out of V R V R . 
This property of Pr turns out to be extremely useful 
in infinite-coordination lattices. In this limit, the con- 
tribution to the average value on |0) of terms in which 
more than two fermionic lines come out of "P R "P R can be 
shown to vanish 14 , which simplifies considerably all cal- 
culations. For instance, the average value on |T]) of any 
local operator Or becomes 

(^V^OrVI^^^vIOrV^), (7) 

which also implies, taking Or = 1, that the variational 
wavefunction ([1]) is normalized. In addition, the average 
value of the inter-site single-particle density matrix turns 
out to be 

(01^4,^,^10) 

= (<l>\-Pk^ a V u V nl c Wt0 V w \<P) 

= ]T ^fz(R) ai Z(R% 5 (0| c Ri7 c n , s |0) 

7(5 

+ E \J Z(TL) ai A(R')^ (5 (0| c R 7 4 (J |0) 

7(5 

+ E \f^UZ(R% s (0| c Rn c R , s |0) 

7<5 

+ £ ^1 A(R) aT A(R')J a (01 c Ri7 4 V |0), (8) 

7<5 

(01^4,04,^10) 
= {<t>\vl^ a v R vl 4,^ R ,|0) 

= £ \/z(RUZ(R<) P s (0| 4, 7 4,5 10) 
+ E y / ^(R) Q7 A(R')^ (0| 4, 7 c R , j5 |0) 

7<5 

+ E \/a(R) Q7 Z(R') i 35 (0| c Ri7 4, )(5 10) 

7<5 

+ £ v/A(R)a 7 A(R0/j5 (0| c R>7 c RV |0), (9) 

7<5 



where the matrices Z and A are determined by inverting 
the following set of equations 

(01^4,^^10) ( 10 ) 

= E \J Z ( K )^ (01 4, 7 C R,/3 10) 

7 

+ E V A ( R ^ (0l c R,7 C R,/3l0)' ( H ) 

7 

(0I^r4,^r4,^I0> (12) 

= e V Z ( R )^ (^ 4,7 4,/j w 

7 

+ E V A ( R )^ (0I C R,7 4, /3 I0) I - (13) 

7 

Naively speaking, it is as if, when calculating the inter- 
site density matrix, a fermionic operator transforms ef- 
fectively into 

4,c = E \I z &)°p 4./3 + V A ( R w c R,^ (1 4 ) 

namely that a particle turns into a particle or a hole 
with probabilities Z and A, respectively. Although all 
the above expressions are strictly valid only in infinite- 
coordination lattices, it is quite common to use the same 
formulas also to evaluate average values on the Gutzwiller 
wavefunction in finite-coordination lattices. This ap- 
proximation is refereed to as the Gutzwiller approxima- 
tioi>Li2 ' 18 i 19 ' 20 , and is known to be equivalent to the sad- 
dle point solution within the slave-boson technique^ 

We conclude by noting that the constraint ([5]) turns 
out to be useful also when the variational wavefunction 
([I]) is applied to Anderson impurity models. In this case 
the operator V acts only on the electronic configurations 
|r) of the impurity, namely 

p = E Arr< |r)(r'|. 

IT' 

If we impose 

(0|ptp ^ = l5 (15) 

(0| V*V C imp \4>) = (0|C mip |0), (16) 

where Ci mp is the single-particle density matrix of the 
impurity, then, for any operator of the conduction bath, 
Obath, and because of (HJ) and (O, the following result 
holds 

(<f>\V 1 O bath V |0) = (0| O bath |0). (17) 

A. Some formal definitions 

In order to perform actual calculations, it is convenient 
to introduce some notations. We define a matrix F a with 
elements 

(F Q ) riF2 = (r i; R| CR jr 2 ;R), 
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as well as its hermitean conjugate, P^, where we assumed 
that the definition of the local configurations is the same 
for all sites, ft follows that 

F a F +F p F a =O, 

where / is the identity. Next, we introduce the uncor- 
rected occupation-probability matrix, P (R), with ele- 
ments 

(P (R)) ri r 2 = (0|ri;R) (T 2 ;R|0), (18) 



that satisfies 



1 = Tr(P (R)), 



Tr (P o (R)FlF 
Tr(p (R)ptpt 



I C R,a C R,/3 



Analogously, the variational parameters that define P R , 
A(R)r 1 r 2 : are interpreted as elements of a matrix A(R). 
With these definitions, Eqs. (UJ) and ([5]) become 

M'Pk'PM = Tr (P (R) A(R)t A(R)) 

= 1, 

= Tr(p (R)A(R) t A(R)FtF /3 



I C R,a c R,/3 



Tr(p (R)A(R)tA(R)F]Fj) 



! c Ra c R,/3 



that suggests to introduce a variational occupation- 
probability matrix P(R) = P (R) A(R) 1, A(R) with ma- 
trix elements 



(P(R)) rir2 = J] (^o(R)) rir3 A(R){, 3r4 A(R) r4lv 

(19) 



r 3 r 4 



that must satisfy 

Tr(P(R)) = 1, 



Tr 







F a ] = (6\cL„c 



R,q ^R,^ 



(20) 
(21) 



(P(K)Fl 

Tr(p(R)ptpt) = (0|4 Q 4 !/3 |0),, (22) 

Eqs. P0]). (f2l~j) and replace the constraints (g]) and 
([5]). With these definitions, the matrices Z and A, see 
Eqs. (fTTjl and (fT3|) . are obtained by solving 



Tr(p (R)A(R)tp t tA(R)P /3 
= ^(R)o 7 Tr (Po(R) P 7 f 



53 JA(R) Q7 Tr (P (R) P 7 PJ , (23) 



Tr(Po(R)A(R)tptA(R)Pj 
= ^ ^Z(R) Q7 Tr (P (R) Pt Pj) 



53 ^(R^Tr (P (R) P 7 Pj) . (24) 



The above equations simplify if one uses the natural 
basis, namely the single-particle basis that diagonalizes 
the density-matrix, 



n(R) a 5 a p, 



In this case 



Tr P (R)A(R)tptA(R)P /3 

= — * ^ — , (25) 

Tr(p (R)A(R)tptA(R)P+) 

A(R) Q/3 - — * : '-■ (26) 

1 - n(R)/3 

Moreover, if one constructs the states |T; R) so that 
Po(R) is diagonal 



(P (R)) rr , =5rr'Po(R;r), 



then 



(P(R)) ri r 2 = PoCRTr) 53 A(R)^ ir3 A(R) r , r3 , (27) 



B. The variational problem 

We are now in position to settle up the variational 
problem. We consider a generic tight-binding Hamilto- 
nian 

n - - V V t°- (i c f c 

RR' a(3 

+ 5353P(R)rr'|r;R)(r';R|, (28) 

R IT' 

where a and /3 stem for spin, orbital and lattice site in 
the chosen unit cell, and the hermitean matrix P(R) with 
elements P(R)rr' may be also unit-cell dependent. The 
average value of this Hamiltonian on the Gutzwiller wave- 
function ([1]) in the limit of infinite coordination lattices 
or, in finite coordination ones, within the Gutzwiller ap- 
proximation, is 



En 



E E 

RR' 

RR' q/37(5 



Z(R) Q7 Z(R% S ■ 



Z(R) Q7 A(R)^< 



"R,7 L R',<5 



I C R, 7 C R',5 



5 



A(R) Q7 Z(R') 



06' 



A(R) Q7 A(R') 



08 



'R,7 C R',<5 



-R,7 L R',<5 



53 Tr(p (R) A(R) f E(R) A(R) 

R 



(29) 



The last term depends only on the local properties of the 
uncorrelated wavefunction \ <f>) , specifically on the occupa- 
tion probabilities Po(R). Therefore, for any given choice 
of Po(R)i the optimal \<f>) that minimizes the variational 
energy is the ground state of the Hamiltonian 



Hvar — 



E 53 

RR' ap-fS 



Z(R) Q7 Z{K')*p & c^ 7 c r , 



Z(R) Q7 A(R); 5 4 7 4 i5 



A(R) Q7 Z(R')^c Ri7 c r , i , 
A(R) Q7 A(R)^ c R7 c R/j 



EE 

R a/3 



fx(R) o/ j 4., « c r„ 



(l>(R)o0 4,a C R,/3 + 



(30) 



where the parameters /i(R) a /3 and i/(R) a/ g are Lagrange 
multipliers to be determined by imposing that the ground 
state has indeed the chosen Pq(R). The last task is to 
find the values of the variational parameters A(R)rr' a s 
well as of Po(R) for which the variational energy (|29]) is 
minimum. We note that the variational Hamiltonian (|30j) 
that has to be solved may include also inter-site pairing 
terms, which are absent in the original Hamiltonian ([28]) . 

Analogously to other more conventional variational ap- 
proaches, like Hartree-Fock theory, it is common to in- 
terpret the single-particle spectrum of the variational 
Hamiltonian (|30[) as an approximation of the true co- 
herent spectrum of quasi-particles.— 



III. THE MODEL 



single-particle density-matrix, a lot of preliminary results 
can be extracted without even specifying how lattice-sites 
are coupled together. Therefore we start our analysis 
from defining some local properties and later we will con- 
sider a specific lattice model. 



A. The isolated dimer 

The basic unit of the model we are going to investigate 
consists of a dimer with Hamiltonian 

U 2 

Udimer = ~t± 53 \ C \a c 1a + H - c -) + ~ 53 ( n * - 



2 ^ 

i=l 



= H± + Hi 



(31) 



where 1 and 2 refer to the two sites of the dimer and n,-, 
i = 1, 2, is the on-site occupation number. 

It is more convenient to work in the basis of the even 
(bonding) and odd (anti-bonding) combinations, defined 
through 

1 / \ 1 , 

Cecr = (CIct + C2a ) , C oa = — = (Cl CT — C2a) . 



V2 



V2 



and use this basis to built the available electronic con- 
figurations, which we will denote as |n, T), with n that 
refers to the number of electrons. The empty and the 
fourfold occupied dimer states are denoted as |0) and 
1 4), respectively, while the singly-occupied states as 

Mo),o) =c + e(o)CT |0), 
and the states with 3 electrons as 

|3, e(o),a) — , , c\ ^ c\ 10). 

1 ' v /' / e(o)a o(e)\ o(e)l 1 ' 

There are six doubly-occupied configurations. Two are 
spin-singlets with two electrons in the even or in the odd 
orbital, |2,e) and |2,o), respectively. When each orbital 
is singly occupied, the two electrons form either a spin 
triplet, |2, l,^) with S z = —1,0,1, or a spin singlet, 
1 2, 0) . Since we are not going to consider variational so- 
lutions that break spin-SU (2) symmetry, it is convenient 
to define the projector operators 

\l,e(o))(l,e(o)\ = 53 \l,e(o),a)(l,e(o),a\, 

a 

|3,e( ))(3,e(o)| - 53 |3, e(o), a) (3, e(o), a\, 



Let us now apply the variational wavefunction to spe- 
cific models that are inspired by the valence-bond crys- 
tal example we mentioned in the introduction, and where 
the off-diagonal elements of the operator Vr as well as 
its non-hermitean character do play an important role. 
Since the operator Pr is built out of purely local prop- 
erties, namely the available on-site electronic configura- 
tions plus a variational guess for the uncorrelated on-site 



|2,1)(2,1| = 



s„=-i 



|2,1,5,)<2,1,5,|. 



The isolated-dimer ground state in the subspace with 
two electrons is 

. . cos 9 / + + \ , . 

I*) = -jf ( c lT c 2i + 4 T c uJ 1°) 



6 



sin 6 

"71" 



( c lT c U + c 2 



2T C 2| 



= —= (cos + sin 0) |2,e) 
v2 

t= (cos 6* — sin0) 1 2 , o) , 

v2 

where tan 20 = it±/U and has energy 



Alternatively, one can consider a paramagnetic |</>) that 
has built in the tendency of each dimer to lock into a 
spin-singlet. This can be accomplished in two ways that 
do not exclude each other. The first is to assume an 
uncorrelated wavefunction with a huge splitting between 
even and odd orbitals, namely with 



£< 



I C L C oa 



E 



u fu 




At 2 , 



(32) 



|\&) can be always rewritten in the form of a Gutzwiller 
wavefunction. First of all, we needs to choose an uncor- 
related wavefunction \<j>). A natural choice might be the 
ground state at U = 0, namely |2,e). Indeed |\E') can be 
written as 



\*)=V\2,e), 



where 



V = |*)(2,e| = (cos (9 + sin (9) |2,e)(2,e| 
v2 



V2 



(cos0-sin0) |2,o)(2,e|. 



(33) 



and obviously satisfies both ^ and ([S]). 

Another possibility, that we are also going to con- 
sider in what follows, is to use an uncorrelated wavefunc- 
tion that corresponds to a dimer in which the two sites 
are only coupled by an intersite singlet-Cooper pairing, 
namely with (c}|Cy) = (c^c^) ^ 0. In this case 

W = 5(|0> + |2 )e >-|2, O >-|4>), 

and, once again, the true ground state can be written as 

|*> = |*> (0| \0)=V\<j>). 

Already at this stage one can appreciate how important 
is the role of the off-diagonal elements in V, especially 
for large U/t±. 



This implies that, among the doubly-occupied configura- 
tions of each dimer, mainly the spin-singlet |2, e) survives 
in the uncorrelated wavefunction. The latter can then be 
turned into the isolated dimer configuration by an ap- 
propriate Gutzwiller operator V, as shown before. The 
other possibility is to include Cooper pairing correlations 
in the singlet channel 



±sc 



I C 1T C 2| 



L 2T C 11 



In this case, the isolated dimer can be recovered by as- 
suming a very strong pairing A$c — 1 and suppressing, 
through V, configurations with none or two singlet-pairs. 
Note that both n e — n a and A$c do not appear by a 
mean- field decoupling of TLu , so that a variational wave- 
function with such correlations built in can not be sta- 
bilized within Hartree-Fock theory. Here the role of V 
becomes crucial. 

Therefore, let us assume for \<j>) a BCS- wavefunction 
defined such that 



cLc 



lazier 



c L C 2a 



C 1T C 2| 



"2(T C 2(T 



C L C la 



n 
Z 4' 
5 

Z 2' 



(34) 
(35) 
(36) 



with real A$c- In the even/odd basis this translates into 

(37) 
(38) 

t uv_ /M JJ uv,^£, (39) 



y 

2 ' 



B. The non-isolated dimer: variational density 
matrix 

When the dimer is coupled to the rest of the system, in 
order to built the operator V we need to specify an un- 
correlated local single-particle density matrix based on a 
variational guess of the uncorrelated wavefunction \<p). A 
simple guess would be a magnetic wavefunction in which 
the two sites of each dimer have opposite magnetization. 
This choice is also the only one admitted by an Hartree- 
Fock decomposition of the interaction term TLu. How- 
ever, a magnetic wavefunction is not the most suitable 
choice to reproduce the limit of isolated dimers, which is 
a collection of singlets. 



where n e + n a = n. As previously mentioned, the cal- 
culations simplify considerably in the natural basis, that 
is derived in the Appendix for this particular choice of 
density matrices. 

As a particular application, we assume hereafter that 
the model is half-filled, namely n e + tiq = 2. The density 
matrix of the operators in the natural basis, d e ^ a and 

d\ r ^_ is, by Eqs. (|A~3)) and (|A~4)) . 



e(o)<7 



where 



d , \ d f x 

e(o)<7 e[o)a 



(40) 
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The two angles 9 e and 8 , that are defined by Eq. (|A.2[) 
and identify the unitary transformation from the original 
to the natural basis, are given by 8 e = 9 and 9 a = — n/2, 
where 



tan26» 



±sc 



(41) 



We note that, for q — > 1/2, the uncorrelated wavefunc- 
tion describes an insulator where charge fluctuations are 
completely suppressed since each natural orbital is fully 
occupied. It is obvious that, if our choice of the varia- 
tional wavefunction is correct, then the optimal uncor- 
related wavefunction must asymptotically acquire 
q = 1/2 for U -> oo. 

The expression in the natural basis of the hopping, 
Eq. (|A.6[) , and interaction, Eq. (|A.7[) . operators can be 
derived through (|A.5|) and have a relatively simple ex- 
pression at half-filling: 



H± = -t± 2cos20(|4)(4| - |0)(0|J (42) 
cos 20 (|3)(3| ~ |T)(T|) - sin 26 (|T)(3| + H.c 
V2 sin 26 (JO) (2, +| + |4> (2, +| + H.c 



U 



|4)(4|-|0)(4|-|4)(0| 



TLu — 



+-(|1)<1| + |3)(3| 
where we have defined 



+[/ (|2, +) (2, +| + |2, -) (2, - 1 + |2, 0) (2, 0| 
U 



(43) 



|T(3)><T(3)| = |l(3),e)(l(3) L e|+|l(3), )(l(3), |, 
|1)(3| - |l,e)(3,e| + |l, )<3, |, 

|2,±) = -L(|2,e>±|2, >), 

and denoted the local configurations in the natural basis 
as |n, r) to distinguish them from the analogous ones in 
the original representation. 



C. The Gutzwiller operator V 

The most general Gutzwiller operator V should include 
at least all the projectors |n,r)(n,r| as well as all the 
off-diagonal operators \n, T) (n', T'\ that appear in the lo- 
cal Hamiltonian, Eqs. (|42f and (|43|) . As we mentioned 
before, our expectation is that the uncorrelated wave- 
function which better connects to the large-C/ Mott in- 
sulator should be identified by q — > 1/2, in which locally 
only the configurations |3) and |4) are occupied with non- 
negligible probability. This suggests that V must include 
at least those off-diagonal operators that would turn |4) 



into the isolated dimer ground state, namely |0)(4| and 
|2, +)(4|. The latter forces to include also |1)(3|, as we 
shall see. 

Therefore we assume for V the following variational 
ansatz: 

V = A„ r |n,r)(n,r| + A 13 |T)(3| 

+Ao 4 |0)(4|+A 2+4 |2,+)(4|, (44) 
with real A's. We define 

P(n,T) = \ 2 nr P (n,r), (45) 

for all n ^ 3, 4, while, for n = 3, 4, 

P(3) = (A§ + A? 3 ) P (3), (46) 
P(4) = {Xl + \ 2 0i + \l +i ) Po(4). (47) 

Then the conditions Eqs. (|4|) and (O read 

J2 P(n, r) = l, (48) 

nP(n,T) = 2 + Aq, (49) 



nr 



A13 Ax VPo(3)P (l) = 



V2A 2+4 A 2+ VPo(4)Po(2, 



(50) 



Here Po(n, T) are the occupation probabilities in the nat- 
ural basis of the uncorrelated wavefunction. Specifically 



■/u(''.n ,/,-,.! ( ^ + qj (- -q 



where g^v is the degeneracy of the configuration. 
Eq. ([50]) guarantees that the anomalous averages 

vanish in the^ natural basis, and explains why we have 
included |1)(3| in (|44|) . It is convenient to rewrite 



; P(4) 
Po(4) 



Ui for i = 4, 04, 2 + 4, 



A, = a — Ui for i = 3, 13, 

V P o( 3 ) 

where it§ + = 1, which can be satisfied by choosing 
U3 = cosip and U13 = sin-0, and U4+Uq 4 -|-u 2+ 4 = 1- The 
latter parameters can be expressed by means of another 
unit vector v = (vi, 1)2,1)3), through 



vi 

V3 



y= («4 + M04) 



cos26> 
sin 20 



(u 4 - m 04 ) - sin26 l M2+4, (51) 



^ (u 4 -uo4)+cos20u 2+ . 
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In terms of all the variational parameters, the P(n,T)'s, 
9, q, ip and v, the average values per dimcr of the in- 
teraction, H-u, and intra-dimer hopping, H±, are readily 
found to be 

Eu = (4\V*HuV\<t>) = ^(P(3) + P(l) 

+ j P(P) + U (p(2, +) + P(2, -) + P(2, 0)) 
+u(v 2 2 +v 2 3 )p(4), (52) 

E± = {cj)\V^n±V\4>) =-2t±S m , (53) 

where the actual correlated values of the hybridization 
and of the anomalous average are 

2,5* = {<j>\V (n e ~n )V\dp) 

= 4 vi v 2 P(4) + cos (29 + 2V>) P(3) 

-2cos26»P(0) - cos26»P(l), (54) 

2A» = 

= 4«iu 3 P(4) + sin(20 + 2V>) P(3) 

-2sin20P(O) - sin20P(l), (55) 

(56) 

We note that <5* and A* are mutually exclusive, namely 
the choice of parameters that maximizes one of the two, 
makes the other vanishing. 

Upon the action of V, the single fcrmion operators 
in the Nambu spinor representation transform effectively 
into 

<p\ ( d e(o) T \ v ^ VZTA e"^ ( d ?< o)T ) , (57) 
V d l(o)i J V d \o)i J 

where Tj, i = 1,2,3, are the Pauli matrices that act on 
the Nambu spinor components, 



tan 3 = —= . 



(58) 



and, finally, 



^dt (o)<T Vd e{o)a \ct>) 



y/P(0)P(l) + -y/P(l)P(2,+) 



+ 2V / P(l)P(2,-)+2V / P(l)P(2,0) 



V3 



+ ^-VP(l)P(2,l) + cos^ ^-^P(3)P(2,1) 



+ cos $ - y/P(S) P(2, 0) + cos ip — y/P{3) P(2, 



+ cos^-VP(3)P(2, 



+ —j= (v\ cos ip + V2 cos (29 + ip) 
v2 v 



+ v 3 sin (20 + if>) V / P(4)P(3) 



(59) 



iptrft u Pd f , v, 

I e(o)T e(o)| 



2- q 



1 



l-4g 2 



- sin^\/P(3)P(2,-) 



1 



^y/P(3)P{2,+) 



sin^ v / P(3)P(2,l) 

1 



+ - smVVP(3)P(2,0) 

+ f«i sin-0 — V2 sin (20 + ip) 
v2 ^ 

+i-3 cos (29 + VO ) VP(4) P(3) 



(60) 



with real V^Z and a/A- Therefore, if the dimers are cou- 
pled one to another by the single particle hopping term 

T= E E 'mr V^'w r (61) 

Rt^R' i,j=e,o 



where 



^R,i — ( c R,i|' c R,i|)' 



and $ its hermitean conjugate, the uncorrelated wave 
function \<p) minimizes the effective hopping 



T var = (Z + A) E E *RR' *R,i T3< 



-2i(3t 2 



R^R' i,j—e,o 



(62) 

under the condition that the local density matrix satisfies 
Eqs. (j3"7} - (j3"!)|) . One can readily show that this amounts to 
find the ground state \<f>) of the variational Hamiltonian 

n var = (z + A) E E ^r^r^^v, 

R^R' i,j=e,o 



R 

E W T l *R,e " *R,o n *R,„) > (63) 



with and /^i such that 
(4>\ H var |0) + 4g// 3 cos(20 + 2(3) + 4g^i sin(20 + 2/3), 



is maximum. 



Before we consider specific lattice models, it is worth 
re-deriving within this variational scheme the isolated- 
dimer ground-state energy (|3"2")l at half-filling. For that 
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purpose, we take all P(n, T) zero but P(4) 
variational energy is simply 



1. The 



E,, 



Err + E I 



U (v 



4i_|_ Vi v 2 - 



The minimum under the constraint v • v = 1 is obtained 
for u 3 = and exactly reproduces (|52")) . We note that 
the minimum energy is independent on 9, namely there 
exists a continuous family of variational solutions with 
equal energy parametrized by 9. However, in spite of the 
fact that the uncorrelated wavefunction may describe a 
superconductor, the actual value of the anomalous aver- 
age A* = 0. 



IV. A LATTICE MODEL OF DIMERS 

As a particular application, let us consider the follow- 
ing lattice model 



n 



53 53 53 tRR ' c ki CT c R', lCT + H c - 



RR' i=l a 

2 



+?EE(»h,«-i) 



R 8=1 



"*JL 53 53 4l/R,2, + H c - 



<RR'> <t 



53 ( £k - *-0 C k,e<7 C k,e<7 + ( £ k + tl 



c k,ocr C k,oCT 



ktr 



u 



(64) 



R i=l 



where 



nR,i = J2<? c R,i - c R,i<7 and £ k is tlie band dis- 
persion induced by £rr', with half-bandwidth D. The 
Hamiltonian represents two Hubbard models cou- 
pled by a single-particle hopping t±, each model being 
defined on a lattice with coordination number z. As we 
mentioned, the variational results that we have so far 
derived are rigorous strictly speaking only if z — > oo, 
although, in the spirit of the Gutzwiller approximation, 
they can be used for generic z as well. 

If U > D,t±, ([64]) describes at half-filling a Mott insu- 
lator which may be magnetic at t± <C D, but is certainly 
non-magnetic at t± D, where the ground state reduces 
essentially to a collection of singlets. For instance, in the 
case of a Bcthe lattice with nearest neighbor hopping, 
the transition is at t± = £)/v°, value that is going to 
decrease if frustration is included. If U is small and the 
Fermi surface is not nested, then the model is metallic 
for tj_ < D and is a band insulator otherwise. In fact, 
in the absence of nesting there is generically a finite win- 
dow of t± values in which, upon increasing U, the model 
undergoes a transition from a paramagnetic metal into a 
non-magnetic Mott insulator, and this is just the case we 
are going to consider in what follows. The same model 
have been recently studied by Fuhrmann, Heilmann and 



Monien using DMFT 23 and by Kancharla and Okamotc 24 
using DMFT and cluster DMFT, respectively, that gives 
us the opportunity to directly check the accuracy of our 
wavefunction. 

The variational Hamiltonian (|63|) of the model (|64|) has 
a very simple expression, 



= 53 *k,e [( £ k* - Ms) T 3 - Tl 



where 



"*k,o [( £ k* + M3) t 3 + m Ti 



£k* = {Z + A) e k . 



k,o> 



(65) 



The variational single-particle spectrum has the conven- 
tional BCS form with eigenvalues 



\J (ek* - M3) 2 + Mi> E ok = \J (ek* + A*3) 2 + 



Mi, 



hence, for any fix ^ 0, has a gap equal to On the 

contrary, when /ii = 0, the spectrum is gapless for | /X3 1 < 
D, otherwise is gaped. The Lagrange multipliers m and 
\xi are obtaining by maximizing 



Ehop — — 53 (-^ek + ^ok) 
k 

+ 4g^ 3 cos(26» + 2/3) + 



1 Hi sin(2(9 + 2/3).(66) 



In terms of (|66p . (|52| and ([53]) the variational energy is 



E r 



Ehop + Ejj + E± , 



(67) 



and depends on eight independent variational parame- 
ters. 

We have solved numerically the variational problem at 
fixed t±/D = 0.5 as function of U/D. To simplify cal- 
culations, we have assumed for the band dispersion ek 
either a flat or a semi-circular density of states, although 
both would give rise to nesting that could stabilize mag- 
netic phases, which we do not take into account. How- 
ever, from the point of view of the paramagnetic-metal 
to paramagnetic-insulator transition, this choice is not 
influential. 

We find that the variational solution displays a first 
order phase transition at U c ~ 2.05 D for a flat density 
of states, as shown by the behavior of the variational en- 
ergy in Fig. [2j This result agrees almost quantitatively 
with the DMFT calculation 2 - obtained with a semicir- 
cular density of states, that also predicts a first order 
transition with a coexistence region between U ~ 1.5 D 
and 1.8 D at the same value of t± = 0.5 D. We note that 
the energy is everywhere finite and vanishes like 1/U for 
large U, see the inset of Fig. [5] The asymptotic behavior 
UE var /D 2 ~ —7/6 is compatible with second order per- 
turbation theory in t and t± using as zeroth-order state 
a collection of dimers, as explained below. In Fig. [3] we 
show the behavior across the transition of the three con- 
tribution to the energy, namely Ejj, E±_ and Ehop- We 
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FIG. 2: Variational energy E vaT in units of D and for a flat 
non-interacting density of states, as function of U/D for t± = 
D/2. At U c — 2.05 D a first order transition occurs. The 
inset shows the asymptotic value of UE var /D 2 . 
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FIG. 3: The different contributions to the variational energy, 
Ejj, E± and Ehop- 



find that the transition is accompanied by an energy loss 
in Ehop, but a gain in both Ej_ and Ejj- 

In order to characterize physically the two phases, in 
Fig. 0] we plot the values of [i\ and /i3 across the tran- 
sition. Since /ii = 0, within our numerical accuracy, 
and I//3I < D, the phase at U < U c is gapless hence 
metallic, see the behavior of the density of states (DOS) 
drawn in Fig. [5] On the contrary, on the U > U c side 



FIG. 4: The behavior of the parameters /ii and fj,3 in units of 
D as function of U/D, see Eq. (I65|l . 



of the transition, \x\ 7^ 0, that implies a finite gap in the 
single-particle variational spectrum, see Fig. [5] In the 
gaped phase at U > U c the spectrum looks like the one 
of a Peierls insulator with a very large hybridization gap, 
not consistent with the bare value of t^. In reality this 
gap is, more properly, the Mott-Hubbard gap. Indeed 
the DOS has weight both below and above the chemical 
potential, suggestive of asymmetric Mott-Hubbard side- 
bands. Moreover, as we are going to discuss below, the 
actual difference between the occupations of the bond- 
ing and anti-bonding bands, which we denoted as 25* in 
Eq. (f54|) . decreases with U, unlike the single-particle gap, 
see Fig. [7J This behavior is reminiscent of what has been 
found by Biermann et alr& as an attempt to understand 
the physics of VO2. 

The other quantities that identify the variational spec- 
trum are Z and A, shown in Fig. [6l We see that Z is 
decreasing with U but reaches a finite value Z — 1/4 for 
U — » 00. On the contrary, A = for U < U c , while 
A 7^ for U > U c and increases monotonically to reach 
asymptotically the same value 1/4 for large U . 

Further insights can be gained by the average values of 
the intra-dimer hopping and pairing, Eqs. (|54|) and (|5*5|) . 
drawn in Fig. [7l As we mentioned, the intra-dimer hy- 
bridization 25* is monotonically decreasing with U, apart 
from the jump at the first order transition. More inter- 
estingly, around the transition the variational solution 
is characterized by a sizeable BCS order parameter A*, 
which does not follow the behavior of the BCS coupling 
Hi present in the variational Hamiltonian. Indeed, while 
fix is zero within our numerical accuracy for U < U C1 
yet a non negligible A* develops just before the transi- 
tion, see Fig. [7] Moreover, although [i\ starts already 
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FIG. 5: The variational single-particle spectrum for the even, 
i.e. bonding, band, solid lines, and odd, i.e. anti-bonding, 
one, dashed lines, across the transition for a non-interacting 
semi-circular density of states. 



large for U > U c and increases monotonically with ?7, 
see Fig. HJ the actual order parameter A* is apprecia- 
ble only near the transition and fastly decreases with 
U to very tiny values. Also interesting is that, besides 
the intra-dimer superconducting order parameter, also an 
inter-dimer one arises, that can be for instance defined 
through 



A 



1 i = l R/ 

4 



0.25 



§eo88 o§ooo< 
A 




_l I I L 



9QP§QQOOQO( 
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U/D 
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FIG. 6: The behavior of the parameters Z and A. 



0.6 - 



0.4 - 




0.2 - 



FIG. 7: The average values of the intra-dimer hopping, 8*, 
and pairing, A*. 



4e k=0 



4V 



R' 



where V is the number of sites. We find that A* has 
actually the opposite sign of A*, and both closely follow 
each other, rapidly decreasing with U, see Fig. |U Even 
though A* and A* are everywhere finite for U > U c , sug- 
gestive of a superconducting phase that survives up to 
very large U, we believe that superconductivity, it it oc- 
curs at all, may appear only very close to the transition, 
where the value of the order parameter is larger. Indeed, 
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FIG. 8: The intra-dimer, circles, and, with reversed sign, the 
inter-dimer, triangles, superconducting order parameters. 



the optimal solution with finite A* and A* and another 
solution in which both are forced to be zero are practi- 
cally degenerate within our numerical accuracy for large 
U. Moreover, since the variational values of the order 
parameter are extremely small but close to U c , the in- 
clusion of quantum fluctuations, for instance in the form 
of a Jastrow factor as in Eq. @, would likely suppress 
superconductivity leading to a bona fide insulating wave- 
function. 

Unfortunately, we can not compare this result with the 
DMFT analyses of Refs. HH1 where superconductivity 
has not been looked for. 26 



A. Large U limit 

In order to appreciate qualities and also single out 
defects of the variational wavefunction, it is worth dis- 
cussing the large-[/ solution. To leading order in 1/U, 
one can assume all P(n, T) = but P(3) and P(4). The 
two constraints Eqs. (|48p and (|49[) can be solved by defin- 
ing 

P(3) = 2 - Aq = Ad, P(4) = 1 - Ad, 

with d <C 1, namely q — > 1/2, while Eq. (f5T))) is already 
satisfied since P(l) = P(2, +) = 0. Moreover, one readily 
recognizes that the variational solution asymptotically 
tends to acquire ^ — /3^7r/4, 0^0 and v% ^> 1)2, V3, 
which is indeed what we find by numerical minimization. 
It then follows that 



•/z 



This implies that /Z3 
by maximizing 



0, hence that fi\ is determined 



E, 



hop 



-Y 

k 



k* 



• Mi 



■4g/ii, 



which leads to pt,\ = y^e^/Ad and 

Ehop = — 4 \J £2 d, 



(68) 



where 



£2 



4V 



At leading order, the variational energy per dimer is 
therefore 



E r 



= -4 V£2 d - 4 t_L V! v 2 + U v\ + 2 U d, 



that is minimized by d — e2/U 2 , v± ~ 1 and t>2 — 2t±/U, 
and takes the value 



2V ^ t/ 

k 



44 

17 



-E 

2V ^ 

RR' 



*RR' ^R'R At 



U 



_J_ 

u ' 



(69) 



In the case of a flat density of states 



k 



£kj 



and with t± — 0.5 73 we recover the numerical result 
E var ~ -7/6 7? 2 /C/, see Fig. HJ 

We note that, in spite of the hybridization <S» ~ t±/U 
being small, the single-particle gap of the variational 
spectrum 2/ii ~ U is large, as one should expect in a 
Mott insulator. 

Coming back to the large-?/ value of the variational 
energy (|69j) , one can readily realize that it coincides with 
the second order correction in iRR< to the energy of the 
state 



I*) = 



n 

R 



i(4u, 



T c R.2i 



C R,2T C R,ll) 



10), 



which is just a collection of singlets. In other words, in 
spite of being non-magnetic, our variational wavefunc- 
tion is able to reproduce the correct super-exchange be- 
tween dimers. This is a remarkable property that ac- 
tually derives from the square-root dependence upon d 
of E hop , see Eq. (|55| . If we considered a more conven- 
tional Gutzwiller operator V commuting with the single- 
particle density matrix, that amounts to further impose 
2q = y/ 5% + A 2 , we would find Eh op oc d, implying a 
transition into an unrealistic insulator with d = above 
a critical U. The obvious defect of the wavefunction is 
that, since it emphasizes strongly the role of individual 
dimers, the hopping among dimers, although finite for 
any U, is under-estimated with respect to the intra-dimer 
one. Therefore we do not expect the wavefunction to be 
particularly accurate for small tj_/D. 
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FIG. 9: The non-interacting density of states of the lattice 
of dimers. The bonding and anti-bonding state of each dimer 
give rise to two bands that overlap, leading to a metallic phase 
in the absence of interaction. 



teraction, to a band insulator. When interaction is taken 
into account, in the form of an on-site repulsion U, one 
expects, above a critical U, a transition from the metal 
into a Mott insulator. If magnetism is prevented, for 
instance by a sufficiently large splitting between bond- 
ing and anti-bonding orbitals and/or by frustration, the 
Mott insulator is non-magnetic. We have shown that 
our wavefunction overcomes the difficulties of Hartree- 
Fock theory and allows to study, albeit variationally, this 
transition. In particular we find that: 

(i) at the variational level the Mott transition is first 
order; 

(ii) the variational spectrum inside the Mott insulator 
looks similar to that of a Peierls insulator with a 
large hybridization gap, namely a large splitting be- 
tween bonding and anti-bonding bands. In reality 
the gap is the Mott-Hubbard gap and the actual 
difference between the occupations of the bonding 
and anti-bonding bands, is small; 



CONCLUSIONS 



(iii) inter-site singlet-superconductivity appears around 
the transition. 



In this work we have proposed an extension of the 
Gutzwiller variational approach to account for correlated 
models which display metal-insulator transitions into 
Mott insulators that escape any simple single-particle 
descriptions, like the Hartree-Fock approximation. The 
wavefunction has still the same form as the conventional 
Gutzwiller wavefunction, 



R 



with R identifying unit cells that may also be non- 
primitive ones, with the novel feature that the opera- 
tor "Pr, is non-hermitean and does not commute with the 
local single-particle density matrix. In essence, this prop- 
erty realizes a variational implementation of a Schrieffer- 
Wolff transformation^, although only restricted to the 
lattice sites within each unit cell. We have shown that, by 
slightly reducing the variational freedom, this wavefunc- 
tion, like the conventional Gutzwiller wavefunction^ al- 
lows for an extension of the Gutzwiller approximation 
to evaluate average values, approximation that becomes 
exact in the limit of infinite coordination lattices. 

As an application, we have considered the Mott tran- 
sition into a Peierls, or valence-bond, insulator, namely 
an insulator that is adiabatically connected to a collec- 
tion of independent dimers. Such an insulator can not 
be described by Hartree-Fock simply because the singlet 
configuration of each dimer is not a Slater determinant. 
Specifically, we have considered the hypothetical situa- 
tion shown in Fig. [91 where the splitting between the 
bonding and antibonding orbitals of each dimer is as- 
sumed not to be sufficient to lead, in the absence of in- 



While (i) and (ii) are presumably true, as they have 
been also found by more rigorous calculation s 23 ! 25 ' 28 , the 
emergence of superconductivity might be an artifact of 
the variational wavefunction^ Nevertheless, the possi- 
ble occurrence of superconductivity is quite suggestive. 
It is known for instance that two-leg Hubbard ladders 
with nearest neighbor hopping display dominant super- 
conducting fluctuations with the same symmetry that we 
find variationally 29 , although at half-filling they always 
describe non-magnetic spin-gaped insulator s 30 ! 31 because 
of nesting. Moreover, the uncorrelated wavefunction 
is quite similar to the wavefunctions used in Refs. l32l 33l 
to simulate t-J ladders. It would be surprising and inter- 
esting if this tendency towards superconductivity turned 
into a true symmetry breaking instability in higher di- 
mensionality, as suggested by our analysis, which we 
think it is worth deserving further investigations. 

Note added: During the completion of this work, we be- 
came aware of a recent extension of slave-boson technique 
whose saddle-point solution closely resembles our varia- 
tional approach. 34 Indeed the two conditions we impose 
on the Gutzwiller operator, Eqs. (g]) and (0, are in one- 
to-one correspondence with the constraints identified in 
Ref. [3J within the slave-boson formalism. 
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APPENDIX: THE NATURAL BASIS 

Let us assume that, in the Nambu-spinor representa- 
tion 



-el 



the uncorrelated wavefunction has the following density 
matrices 

A _( n e /2 A S c/2 \ a _( n /2 -A 5C /2 
° e - { A sc /2 1 - n e /2 )> Uo ~{ -Asc/2 1 - n /2 

(A.l) 

The natural orbitals are obtained by the unitary trans- 
formation 



d e(o) T = COS0 e(o)C e{o)1 +Sme e(o) cl (o)l 



e(o)l 



COS I 



e(o) C e(o)l 



sin I 



e(o) C e(o)T' 



where 



tan2# P = — tan20 o = 



n R - 1 



n - 1 



(A.2) 



and posses a diagonal density matrix with the non- 
vanishing elements given by 



1 



Qe(o) 



(A.3) 



where 



<?e(c) - ^ ^(ne^-l^ + Alc. (A.4) 

In the natural basis we introduce states that have the 
same formal expression as in the original basis but are 
built with d-operators, and denote them as \n,T). The 
transformation rules from these states to the original ones 
are 

|0) = cos6> e cos0 o |O) +sin# e sin# |4) 

+ cos£? e sin^o 1 2, o) + sin 9 e cos0 o |2, e), 
\l,e(o),a) = cos6» o(e) |T, e(o),a) + sin6> o(e ) |3, e(o),cr), 



|2,e(o)) = cos0 e cos0 o |2, e(o)) + cos0 e ( o ) sin^ ( e ) |4) 
-sin0 e(o) cos6> ( e) |0) 
-sin6> e(o) sin6> o(e) |2, o(e)), 
|2,1,S,) = |2,1,S,), (A.5) 
|2,0) = |2,0), 
|3,e(o),cr) = cos6» o(e) |3, e(o)) - sin6» o(e ) |l,e(o)cr), 
1 4) = cos6» e cos6> |4) +sin6» e sin6> |0) 

— cos£? e s'm8 1 2, e) — s'm9 e cos 9 Q |2, o). 

The inverse transformation is obtained by letting 9 e r -\ — > 

-#e(o)- 

The hopping operator in the original representation 
can be written as 



H± = -2t ± c l c 2* + H.c. 

a 

= -2t± C L C ea - C L C e 



-2t i [3,0,17) (3, o, crj - |3, e, a) (3, e, a 

cr 

+ |l,e,cr)(l,e,er| - \l,o,a)(l,o,a\ 
+2|2,e)(2,e| -2|2,o)(2,o 



(A.6) 



while the interaction operator as 

2 



n 



u 



u 



i=l 



- = *L 

2 



>(0| + |4)(4| 



Zl X/ l n > e '°')( rl ' e > 'l + l n )°! Cr )( n )°! 'l 

a n=l,3 

-(|2,e> + |2, >) ((2,e| + (2, f 



-2|2,0)<2,0| 



(A.7) 



Their expression in the natural basis can be obtained by 
the transformation rules (|A.5|> . 
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